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EFFECT OF CONTROLLED SURFACE ROUGHNESS 


ON BOUNDARY- LAYER TRANSITION AND HEAT TRANSFER 
AT MACH NUMBERS OF 4.8 AND 6.0 
By Paul F. Holloway and James R. Sterrett 


SUMMARY 


An investigation has been conducted at a free -stream Mach number of 6.0 to 
determine the effects of controlled three-dimensional surface roughness (spheres) 
on boundary-layer transition and heat transfer. Experimental data are presented 
for a sharp-leading- edge two-dimensional flat-plate model over a local unit 
Reynolds number range of 1.2 x 10^ to 8.2 X 10^ per foot based on a local Mach 
number of 6.0 and 9*8 X 10^ to 13-9 X 10^ per foot based on a local Mach number 
of 4.8. 


The Reynolds number for natural transition has been found to increase mark- 
edly with increasing Mach number above a Mach number of approximately 3-7- It 
was found that surface roughness of height less than the boundary-layer thick- 
ness can delay transition. The critical roughness Reynolds numbers determined 
experimentally in this investigation are larger than those found at lower super- 
sonic Mach numbers. Roughness heights of approximately twice the calculated 
boundary- layer thickness at the roughness position were required to trip the 
boundary layer completely. Calculations of the heat-transfer distributions based 
on simple flat-plate theory are shown to agree reasonably well with the exper- 
imental results . 


INTRODUCTION 


The nature of boundary- layer transition and its effect on associated problem 
areas such as surface heat-transfer rates and aerodynamic characteristics has 
been the subject of many studies, yet there is still much to be learned about the 
phenomena. The data and theories on boundary- layer transition currently avail- 
able for higher supersonic and hypersonic Mach numbers are limited at best and 
contain a great deal of scatter. Indeed, the designer of a high-speed configura- 
tion finds that it is difficult to evaluate the meaning of the available con- 
tradictory experimental and theoretical results. 

One of the most important factors affecting the transition of the boundary 
layer from laminar to turbulent is the condition of the body surface. Surface 
irregularities such as expansion slots envisioned as necessary in winged reentry 



vehicles (for example, the X-20), connecting rivets, or distortions due to high- 
temperature buckling of the skin material may well be prominent factors in deter- 
mining the nature of the local boundary layer. Since it is possible that surface 
irregularities may cause boundary-layer transition, the designer of winged 
reentry configurations must be able to estimate the extent to which these irreg- 
ularities may affect the location and length of transition and, in turn, the con- 
sequential effect on the surface heat-transfer rates. 

Surface condition is also of importance from the experimental point of view 
since roughness may be utilized as a boundary- layer trip in wind-tunnel or free- 
f light tests. A guide is needed for determining the size of roughness required 
at high supersonic and hypersonic Mach numbers to obtain fully developed turbulent 
boundary layers which will duplicate full-scale conditions. 

Considerable study has been devoted to boundary- layer transition and the 
factors affecting transition at subsonic and supersonic speeds. (See, for 
example, refs. 1 to 7- ) In addition, recent work on the effects of surface dis- 
tortions on the heat transfer to a wing at hypersonic speeds has been published 
in reference 8. 

The purpose of this report is to present the results of an experimental 
study of the effects of one type of surface distortion of importance, controlled 
surface roughness (spheres), on boundary-layer transition and surface heat- 
transfer rates at Mach numbers of ^.8 and 6.0. The experimental investigation 
was conducted on flat plates with sharp leading edges in the Langley 20- inch 
Mach 6 tunnel and in the Mach number 6.2 blowdown tunnel of the Langley Research 
Center. The results are compared with the available data, empirical results, 
and theory from the literature for supersonic and low hypersonic Mach numbers. 


SYMBOLS 


Ai, A 2 

Cf 

Cp 

c p,o 


C W 


d 

h 

k 


coefficients of T f equation 

local skin-friction coefficient 

average skin- friction coefficient 

specific heat at outer edge of boundary layer 

specific heat of wall material 
diameter of roughness elements 
heat-transfer coefficient 
vertical height of roughness above plate 
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I 

M 

N Pr 

%t 

4 

r 

% 

Ro 

R o,t 


Rr 


*0,X* 


R x,T ! 


s 

t 

T 

T ' 

u 


length of transition region 
Mach number 
Prandtl number 

Stanton number 
experimental heating rate 
recovery factor 

Reynolds number based on fluid conditions at top of roughness elements 


and height of roughness. 




unit Reynolds number per foot at outer edge of boundary layer. 


H o o 

Ho 


transition Reynolds number. 


PqV 

Ho 


local free- stream Reynolds number based on distance from virtual 

Pq^o^v 
origin, „ 

HO 


local free-stream Reynolds number based on distance from leading 
edge. 


Pq u o x 

Ho 


local free-stream Reynolds number based on distance from roughness 
p u x* 
location, - ° — 

H 0 

Reynolds number based on T 1 conditions and distance from virtual 

. . P u o x v 
origin, 

lateral spacing of roughness 

distance from leading edge to transition (either at end of laminar 
flow or at beginning of turbulent flow) 

temperature 

reference temperature 

velocity component of flow parallel to surface of plate 
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x distance from leading edge 

x k distance from leading edge to roughness position 

x v distance from virtual origin 

x* distance from roughness to instrumentation location 

y vertical coordinate measured from plate surface 

p density 

fi leading-edge thickness 

9 W local vail thickness 

t time 

6 calculated boundary- layer thickness at roughness position based on 

velocity 

a angle of attack 

7 ratio of specific heats 

\x viscosity 

Subscripts: 

1, 2,3 denote relative roughness height 

cr critical 

k conditions at top of roughness elements 

lam laminar conditions 

o local conditions at outer edge of boundary layer 

r recovery 

T 1 based on reference temperature conditions 

trans transitional conditions 

turb turbulent conditions 

v vail 


If 



00 


free- stream conditions 


Primes denote parameters evaluated at reference temperature T'. 


APPARATUS AND TEST METHODS 


Wind Tunnel 

The test program was conducted in the Langley 20- inch Mach 6 tunnel and in 
the Mach number 6.2 blowdown tunnel. The Langley 20-inch Mach 6 tunnel is the 
intermittent type exhausting to the atmosphere through a diffuser augmented by 
an air ejector. Tests were run with tunnel stagnation pressures of 3^5; 440, 
and 515 pounds per square inch absolute with stagnation temperatures of 96 O 0 R 
to 1,020° R. A more detailed description of the tunnel is given in reference 9. 

In order to extend the test Reynolds number range below that obtainable in 
the Langley 20-inch Mach 6 tunnel, additional tests were conducted in the Mach 
number 6.2 blowdown tunnel. The tunnel is also of the intermittent type 
exhausting to a 40, 000-cubic-foot sphere which can be pumped to pressures as low 
as 1 millimeter of mercury absolute. Tests were run with tunnel stagnation pres- 
sures of approximately 65 ^ 165 , and 265 pounds per square inch absolute with 
stagnation temperatures of 840° to 1,020° R. A more detailed description of the 
tunnel is given in reference 10. 


Models 

The models tested consisted of flat plates with sharp leading edges con- 
structed from stainless steel. Each model assembly was 9 inches wide and approx- 
imetely 11 inches long for the tests in the Langley 20-inch Mach 6 tunnel. 

Plate 1 was a continuous plate, a sketch of which is given in figure l(a). The 
remaining models consisted of plate 2 with interchangeable leading edges (one 
for each roughness height). (See fig. 1(b).) The instrumented plates and the 
leading edges were mounted on a support plate as shown in figure 1. The smaller 
size of the Mach number 6.2 blowdown tunnel required that the models be smaller. 
For the tests in this tunnel, plate 2 and the leading edges were cut down so that 

the model assembly was 7^ inches wide and loi inches long. 

For the roughness tests, the leading-edge pieces were interchanged, each 
having a different size roughness mounted 2 inches from the leading edge. Flat 
plate 2 and the support plate completed the assembly. (See fig. 1(b).) The 
spheres were glued into small spherical indentations in the leading-edge plate. 

The location, spacing, height above the plate, and diameter of the spheres are 
given in figure l(b). 

The leading edge of all the models was a 20° wedge that tapered to a cylin- 
drical leading edge with a radius of approximately 0.002 inch or less. Two models 



of each plate were constructed - one being instrumented with 0.050-inch ID pressure 
orifices and the other with 30-gage iron- constant an thermocouples. The instrumen- 
tation was located chordwise along the center line of the plates. The undersur- 
face of each plate instrumented with thermocouples was slotted along the center 
line to a width of 0.6 inch and a surface skin thickness of approximately 
0.020 inch to minimize the lateral heat conduction in the skin. 


Test Methods and Techniques 

For the tests in the Langley 20-inch Mach 6 tunnel, the free- stream Mach 
number outside the boundary layer of the models was varied by changing the angle 
of attack of the plates. This method gave a local surface Mach number M 0 of 
6.0 at an angle of attack of 0° and of 4.8 at an angle of attack of 8 °. The 
resulting unit Reynolds numbers per foot based on conditions outside the boundary 
layer were 5.82 X 10^, 7-02 X 10^ and 8.22 x 10^ for Mo = 6.0 and 9*84 X 10^, 
11.86 X 106^ and 13*92 X 106 for Mo - 4.8 for the tunnel stagnation pressures 
of 365 , 440, and 515 pounds per square inch absolute, respectively. 

The tests were run in the Mach number 6.2 blowdown tunnel with the plates at 
an angle of attack of 0°, a nominal local Mach number of 6.0 being obtained. The 
resulting local unit Reynolds numbers per foot were approximately 4 x 106 ^ 

2.6 X 106, and 1.2 X 106. 

Pressure test s.- Pressure distributions along the center line of the plates 
were obtained in the Langley 20-inch Mach 6 tunnel for Mo = 6.0 and 4.8 with 
free- stream Reynolds numbers per foot of 5.82 x 106^ 7 . 02 X 106^ and 8.22 X I 06 . 
The local static pressures on the plates were measured by connecting the orifices 
to pressure-switching devices which in turn connected the orifices in sequence to 
electrical pressure transducers. The electrical outputs from the transducers were 
recorded on a digital readout recorder. Each pressure- switching device was con- 
nected to two transducers with ranges of 1 and 5 pounds per square inch absolute. 
The accuracy of the transducers is approximately l/2 percent of full-scale 
reading. All pressure tests were run on the same support system as was used for 
the heat-transfer tests. 

Heat-transfer tests .- The aerodynamic heating was determined by the transient 
calorimetry technique by which the rate of heat storage in the model skin is meas- 
ured. The models, originally at room temperature or slightly cooler, were sud- 
denly exposed to the air flow by quick injection from a sheltered position beyond 
the tunnel wall. Injection was accomplished in less than 0.25 second and the 
model remained in the tunnel for only 4 seconds so that the model was in a nearly 
isothermal condition, lateral conduction in the skin being kept to a minimum. 

Optical methods .- During the pressure and heat-transfer tests in the Langley 
20-inch Mach 6 tunnel, shadowgraphs and schlieren photographs were occasionally 
taken to aid in determining the type of boundary layer existing on the plates. 
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DATA REDUCTION 


The electrical outputs from the thermocouples were recorded on a high-speed 
digital readout recorder. The reading from each thermocouple was recorded at 
0 . 025 - second intervals, converted to a binary digital system, and recorded on 
magnetic tape. The temperature- time data were fitted to a second-degree curve by 
the method of least squares, and the time derivative of temperature was computed 
on a card-programed computer. 

The tunnel stagnation temperature range was 8k-0° R to 1,020° R and the wall 
temperature of the plate was approximately 55 O 0 R. Because of the short time 
required for the injection of the model, the plates were considered to have been 
subjected to a step function in the applied heat-transfer coefficient. The thin- 
skin equation used to calculate the local surface heating rate was 

, hT-rj f * 

q = c w p w 9 w (l) 


The measured local heat-transfer coefficient was then calculated by the relation 


h 



( 2 ) 


in which conduction effects are neglected, and where T r is the calculated 
recovery temperature defined as 


T r = T r 


1 + Mo 


2 r(/ - 1 ) 


T w is the measured wall temperature, and Mq is the local Mach number outside 
the boundary layer calculated from the measured pressure distribution. For the 
tests in the Mach number 6.2 blowdown tunnel, the pressure distribution was 
assumed to be the same as that obtained in the Langley 20-inch Mach 6 tunnel for 
a given configuration. The recovery temperature T r was calculated by assuming 
a recovery factor equal to O. 83 O in the laminar region and O .883 in the turbulent 
region. For the transitional region, T r was calculated by assuming a linear 
variation given by 


( 3 ' 


where 

Ax linear distance from beginning of boundary- layer transition 

l length of transition region 

Finally, the Stanton number, based on local conditions outside of the boundary 
layer, was calculated by use of the equation 
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h 


00 


N< 


St = 


p o u o c p,o 


The experimental heat-transfer parameters q, h, and Ng-^ presented in 

this report were determined by reading the slope of the temperature- time curve at 
a time 0.05 second after the model was in position in the tunnel. The maximum 
surface temperature increase at the time for which the parameters were calculated 
was always less than 25° and generally less than 15°. This low temperature 
increase combined with the thin-skin thickness kept the conduction error at a 
minimum. The inaccuracy of the distribution is thought to be less than 
±10 percent. 


REVIEW OF HEAT- TRANSFER EQUATIONS 


There are many methods available for the theoretical calculation of the 
Stanton number. (See, for example, refs. 11 to 16.) In this analysis, the T T 
method of Monaghan (refs. 13 and 14) has been employed. From reference 1 6 , the 
T* equation may be written as 


T ! = ApT w + T 0 


1 - A 1 + a 2 ^21_ — 


(5) 


The Stanton number may be determined from the modified Reynolds analogy 


N 


St,T‘ - 2 


1 fy T ' » 


(Np r ) 


2/5 


(6) 


where the Prandtl number Np r is based on T* conditions, 


Laminar Boundary Layer 

For the case of a laminar boundary layer, the coefficients of equation ( 5 ) 
(see ref. 12 ) become 




Ap = 1 - O.kSd^N-pj,) 1 ^ 

A 2 = (l - A x - 0.273N Pr )(N Pr ) 1 / 2 

By using the Blasius equation for laminar flow, equation ( 6 ) becomes 

0.332 ^C 7 


(7) 


N 


St 


^x^Pr ) 2 / 5 


( 8 ) 
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where Ngt is based on free-stream conditions for direct comparison with data, 
and the conversion parameter from T' reference conditions to free-stream con- 
ditions C' is given by 


C* 


|i’ T 0 
Mo T' 


(9) 


Turbulent Boundary Layer 


For the case of a turbulent boundary layer, the coefficients of equation (5) 
(see ref. 13) become 


( 10 ) 



The Karman-Schoenherr equations were used to determine the local skin-friction 
coefficient as follows. (See ref. 15 for a plot of these parameters and further 
discussion.) 


0.242 

\[ C F, T ' 


" log lo( R x,T , )( C F,T') 


( 11 ) 


where 


Rx,T' 


P’ u o x v 

M' 


( 12 ) 


and C-^ ip t 


C f ,T 1 = 1 

C F,T' 1 + 3-59\/ c F,T’ 

may be converted to free-stream conditions as follows: 


Cf 


C f,T' 

T’/Tq 


(13) 


(14) 


and finally Stanton number is determined by modifying equation (6), basing Ng-j. 


on free-stream conditions so that 

%t 


1 

2 


Cf 

(Npr ) 2/5 


(15) 


where Cf is determined by equation (l4-). The distance x v of equation (12) 
is defined as the distance from the virtual origin. In this paper, the virtual 
origin for the flat plate with undisturbed flow was defined as the point at which 
laminar flow ends, as determined by the surface heat-transfer rates. (See 
ref. 15 for further discussion.) For the flat plate with roughness, the virtual 
origin was found by determining the length of transition l on the smooth plate 
for a given free-stream Reynolds number, assuming fully developed turbulent flow 
to begin at the trip position, and then defining the virtual origin to be located 
at a distance l forward of the trip position. 
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For the lower Reynolds number tests, fully developed turbulent flow was not 
obtained on the smooth flat plate; therefore, the length of transition was not 
known. In order to compare the experimental results with theory for these tests, 
two theoretical distributions were calculated by: 

(1) Assuming the virtual origin to be located at the roughness position 

( R o,v = R o,x*) 

(2) Assuming the Reynolds number based on the virtual origin to be 

R 0 ^ v = Rq^x* + 2.7 X 10^, where 2.7 x 10^ was based on the length of 
transition for the higher Reynolds number tests. 


Transitional Boundary Layer 

A simple semi-empirical method of predicting the heat transfer in the transi- 
tional boundary- layer region has been presented. This method is based on the near 
linear increase in heat transfer which begins at the end of laminar flow and is 
presented only for general comparison purposes. In this method, it is assumed 
that the Stanton number increases linearly from the end of laminar flow to the 
point at which fully developed turbulent flow is first obtained. The difficulty 
with the relation is that the beginning and end of transition of the boundary 
layer must be determined before it can be applied. 


RESULTS AN D DISCUSSION 


Determination of Transition 

As has been discussed by Probstein and Lin in reference 2, many factors are 
known to influence the selection of the so-called "transition point" from laminar 
to turbulent flow in wind-tunnel tests. For instance, the transition process 
occurs over a finite distance covering a significant range of Reynolds numbers; 
and the type of instrumentation and method used to detect transition influences 
the selection of the transition point. 

An accurate method of detecting transition is needed for the proper evalua- 
tion of many wind-tunnel test results. In this report, the location of transi- 
tion has been determined primarily by noting a change in the local surface heating 
rate along the plates. A typical example of this effect is shown in figure 2 
which gives the distribution of the local heating rate along the plate at two 
local Mach numbers and several unit Reynolds numbers. As the laminar boundary 
layer on the forward end of the plate thickens, the local heating rate decreases 
until transition occurs. Once transition of the boundary layer begins, the local 
heating rate increases rapidly. When fully developed turbulent flow is obtained, 
the heating rate peaks and begins to decrease with increasing distance from the 
leading edge. The general shape of the surf ace-heating- rate distribution is, of 
course, similar to that of the usual shear distributions for a flat plate. (See, 
for example, ref. 17 •) The locations defined as the beginning and end of the 
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transition region are sketched on the insert of figure 2(a). The beginning of 
transition (end of laminar flow) was taken as that region where the heating rate 
begins to increase rapidly. The end of transition (beginning of fully developed 
turbulent flow) was taken as that region where the local heating rates begin to 
level off rapidly. An objection to the surf ace-heating- rate method of deter- 
mining transition is the error introduced by lateral thermal conduction in the 
skin. However, in these tests great care has been taken to minimize the thermal- 
conduction in the skin so that the transition data presented herein are thought 
to be accurate. 

Mach number, Reynolds number, surface roughness, leading-edge thickness, 
temperature of model, free- stream turbulence, and so forth, all affect the loca- 
tion of transition. Also, the method of detecting transition influences the 
apparent experimental location of transition. Therefore, as might be expected, 
a compilation of the available transition data shows considerable scatter. These 
data were obtained from references l 8 to 25 and the present investigation and are 
presented in figure 3 with transition Reynolds number plotted as a function of 
Mach number. Reference 25 has indicated that two very important parameters for 
correlating the data from the various installations are free- stream unit Reynolds 
number and leading-edge thickness. However, a further correlation of the data 
in figure 3 is difficult and out of the scope of this paper; rather figure 3 is 
meant to show that transition data of this investigation agree reasonably well 
with those found in the literature and that the general trends of the available 
data indicate that the stability of the boundary layer increases with increasing 
Mach number above Mach numbers of approximately 3*5 to 4. Similar results have 
been found by Potter and Whitfield in reference 2 6 where for a given unit Reynolds 
number and leading-edge thickness, the transition Reynolds number was found to 
increase markedly above a free-stream Mach number of approximately 4.0. The 
transition Reynolds number is shown (in ref. 26) to increase as much as 500 per- 
cent in going from a free-stream Mach number of 4.0 to a free-stream Mach number 
of 8 . 0 . 

To complement the data determined by the heat-transfer method, shadowgraph 
and schlieren photographs (see, for example, refs. 9 and 27 ) were also employed 
to detect the location of transition. A typical shadowgraph for the flat plate 
is shown in figure 4(a). The beginning of a change in the slope of the thick 
white band represents the beginning of transition. The end of transition is 
indicated when the band converges to the apparent plate surface. A typical 
schlieren photograph of the flow over the flat plate is shown in figure 4 (b). 

The beginning of transition is assumed to occur at the point on the schlieren 
photograph where the boundary layer begins to thicken more rapidly. (See 
ref. 18. ) A comparison of the transition data from the shadowgraph and schlieren 
photographs with no roughness with that from heat-transfer measurements (fig. 2 ) 
as given in figure 3 shows reasonable agreement in the location of transition 
by the several methods when it is remembered that the beginning of transition 
is not a point but rather a small region which may shift somewhat with time. 

As was pointed out in reference 2, the work of Coles (ref. 20) shows that 
the Reynolds number for the end of transition differs from the Reynolds number 
for the beginning of transition by a factor of 1.5 for supersonic Mach numbers. 

It is interesting to note (from fig. 3 ) that the data of this investigation fit 
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approximately this empirical relation for a given free- stream Reynolds number 
at local free- stream Mach numbers of 4.8 and 6 . 0 . 


Effects of Roughness on Boundary- Layer Transition 

One of the governing variables in producing transition by three-dimensional 
roughness elements in supersonic flow is the height of the roughness. Figure 5 
presents a sketch typical of the transition positions for various height spheres 
taken from the data of references 5 and 6 for supersonic flow with zero heat 
transfer. The different regions of transition as indicated by the curves for 
various sized roughness are labeled 1 , 2 , and 3 * In zone 1 , the free-stream 
disturbances are predominant in establishing transition, whereas the disturb- 
ances from the roughness play the predominant role in zone 3 and a further 
increase in sphere height has little effect on the location of transition. In 
zone 2, both the free-stream disturbances and the disturbances created by the 
roughness have an effect on the position of transition, and a further increase in 
the height of the roughness will cause the transition location to move forward. 
The data of reference 5 also indicate that the lateral spacing of a single row 
of spheres has little effect on boundary-layer transition provided the spheres 
are not so close together that they act as a two-dimensional roughness element. 
The effect of sphere spacing was not examined in this test program; however, a 
roughness element was always located on the center line of the model, forward of 
the thermocouples (fig. l) in order to minimize any effect of sphere spacing. 

In this paper, a roughness trip is defined as effective (or critical) when 
a further increase in roughness height causes little change in the forward move- 
ment of the location of the end of the boundary- layer transition. 

In figures 6 and 7 , the distributions of the local heating rates along the 
plate for various height roughness are presented for local free-stream Mach num- 
bers of 6.0 and 4 . 8 , respectively. In these figures, the data given by the cir- 
cular symbol represent the heating rates obtained on plate 2 with a sharp leading 
edge and no roughness. The dashed line represents the faired data from figure 2 
for the continuous plate 1 . The discrepancy between the two sets of data is 
thought to result not only from small variations in the leading-edge thickness 
(see, for example, refs. 25 and 26), but also from a small angle -of- attack varia- 
tion (less than 1 / 2 °) between the two assemblies. However, the angle of attack 
for each series of roughness tests was invariant since the data represent a 
group of tests made without changing the mounting. 

Examination of figures 6(a) and 6(c) shows that the smaller diameter rough- 
ness elements (k = 0.0018 foot and 0.0030 foot) were actually found to delay 
transition beyond the transition point found with plate 2 and no roughness. 

Since these data indicated that small roughness heights may delay transition, 
further limited roughness tests were conducted with the model assembly having 
the same leading edge so that possible effects of slight variations in nose 
bluntness would be eliminated. (This work was stimulated by the comments of 
Potter and >/hitfield In ref. 28 in regard to the summary of this study published 
in ref. 29.) Figure 6(f) presents the heat flow rate distributions for the 
model with the same leading edge, and with various height roughness elements at 
a unit Reynolds number similar to that of figure 6(a). A schematic of the model 
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shown in the figure illustrates the new mounting technique, the roughness strips 
being interchangeable and with the same leading edge (thickness diameter of 
0.0025 inch), being used for all data of this figure. Comparison of the data in 
figures 6 (a) and 6 (f) indicates that although some of the variation in the loca- 
tion of transition with small roughness and no roughness was due to the slight 
variations in the leading-edge thickness, under certain conditions transition is 
apparently slightly delayed when the surface roughness is less than the boundary- 
layer thickness. 

The physical reasons for these phenomena are not fully understood; however, 
reference 19 has indicated that the stability of a separated laminar mixing layer 
increases with Mach number more rapidly than the stability of an attached laminar 
boundary layer. This fact has led to the speculation that the laminar separation 
that must exist near the small roughness elements is at least partially respon- 
sible for a delay in the transition at the present Mach number. However, the 
pressure loss at the edge of the boundary layer caused by the roughness elements 
may be a contributing factor. These phenomena appear to warrant further 
i nve s t igat ion . 

For the higher Reynolds number tests at Mq = 6.0 (figs. 6 (a), 6 (b), and 
6 (c)) as the roughness is increased above k = 0.003 foot, transition moves for- 
ward until a further increase in height is ineffective. For example, note that 
in figure 6 (a) the heating- rate distributions for the roughness of height 
k ^ 0.0054 foot are approximately the same, and k = 0 . 005 * 4 - foot is defined as 
the roughness height slightly greater than the critical roughness height for 
this free-stream Reynolds number. 

Similar trends are evident in figure 7 ^ or ^o = 4.8, the smallest size 
roughness (k = 0.0018 foot) now causing a considerably greater forward movement 
of the beginning of transition than was found for M 0 = 6.0, since the local 
Reynolds number is increased. However, the smallest size roughness is still 
below the critical height as previously defined. 

Previous work has established a criterion for determining the critical 
height of three-dimensional roughness elements for subsonic and supersonic 
speeds. (See, for example, refs. 4 and 30.) This criterion states that transi- 
tion will occur when the roughness Reynolds number exceeds a certain value which 
is approximately constant with the Mach number. This critical roughness Reynolds 
number is, by definition, based on the fluid properties at the top of the rough- 
ness element and the height of the roughness element. The value of the critical 
roughness Reynolds number was not found to be affected appreciably by moderate 
surface cooling. (See ref. 30*) 

The variation of roughness Reynolds number R& with various assumed rough- 
ness heights is presented in figures 8 and 9 ^ or local free-stream Mach numbers 
of 6.0 and 4.8, respectively. These figures show the value of for both 

adiabatic wall conditions (T w /Tr = 1 . 0 ) and the actual wall temperature 
(Tv/T r = 0.66) of the present tests. The roughness Reynolds numbers presented 
were determined from boundary- layer temperature and velocity profiles calculated 
by the Chapman-Rubesin method for an x-distance of 2 inches from the leading edge 
(the roughness location). A complete discussion of the validity and assumptions 
of this theory is given in reference 31 . The velocity profiles calculated by 
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this method are also given in figures 8 and 9- In addition, these figures show 
that the effect of cooling the boundary layer is to cause a thinning of the 
boundary layer and thereby increase the Rk values for a given height roughness 
provided the roughness height is less than the boundary- layer thickness. 

The critical roughness Reynolds number has been obtained experimentally for 
lover supersonic Mach numbers by several investigators. A summary of the experi- 
mental results from reference 4 (-which included data from other sources, 
i ncludi ng ref. 5) is presented in figure 10 in the form of the variations of 

cr with Mach number. Also plotted in this figure are the results of the 
present investigation where the open symbol indicates that the roughness height 
is slightly less than the critical height and the solid symbol indicates that 
the roughness height is slightly greater than the critical height. This figure 
shows clearly that the critical roughness Reynolds numbers determined experi- 
mentally for Mo = 4.8 and Mq = 6.0 are much larger than those observed in 
reference 4 for lower supersonic velocities. The critical roughness height of 
the present tests is also much larger than that predicted by the semi-empirical 
equation of reference 6, which was determined from experiments conducted at 
lower supersonic speeds. 

It should be mentioned that the definition and determination of the critical 
Reynolds number for the previous data shown in figure 10 from reference 4 are not 
exactly the same as that used in the present report. Also, the data of refer- 
ence 4 were taken under zero heat-transfer conditions. Some possible values 
based on the end of transition from the data of references 2 6 and 32 are also 
shown in figure 10. In reference 32, the end of transition was determined by 
the shadowgraph method which is similar to the present transition detection 
method since both methods detect permanent changes in the boundary- layer pro- 
files. The data in figure 10 show that basing the critical Reynolds number on 
the end of transition gives values which are higher than those given in 
reference 4. 

A modification of the critical roughness correlation parameter R^ has 
been given by Potter and Whitfield in references 2 6 and 33* The results of the 
present investigation agree well with the extrapolated values given in these 
references where it is predicted that for local hypersonic Mach number and walls 
at temperatures corresponding to the conditions of the present investigation, 
the Rk value needed to bring about transition at the roughness elements would 
be approximately 2.4 x 10^. The values of Rk of the present investigation 
were approximately 2 x 10^" to 4 X 10^. 

The experimental variation of the critical roughness Reynolds number with 
free-stream unit Reynolds number for a Mach number of 6.0 is shown in figure 11. 
Included in this figure are both the data taken with one leading edge with rough- 
ness elements mounted on an interchangeable strip and the data taken with inter- 
changeable leading edges. The two sets of data do not coincide; however, this 
condition may be due in part to a different location xk of the roughness ele- 
ments. An increase in free-stream unit Reynolds number by a factor of approxi- 
mately 8 was found to result in an increase in critical roughness Reynolds num- 
ber by a factor of approximately 4. Hence at a Mach number of 6.0, cr has 
been found to be sensitive to the free-stream unit Reynolds number. 


14 



For local free-stream Mach numbers of 4-. 8 and 6.0 , it has been shown herein 
that the roughness parameter k/b must be approximately 2 or greater to move 
the beginning of fully developed turbulent flow to the region of the roughness. 

As is pointed out in reference 26, this requirement of a roughness height of 
twice the boundary- layer thickness may lead to limitations in the use of three- 
dimensional roughness to obtain turbulent flow in wind-tunnel tests because of 
the flow distortions created by the roughness extending outside the boundary 
layer. Of course, the requirement of k/5 values of 2 or greater is partially 
responsible for the large variation in Rj^ cr "with varying free-stream Reynolds 
number. 

Further work is warranted to give a better insight into the importance of 
the defining parameters in boundary- layer transition work. It is interesting to 
note that the recent theoretical investigation of Lees and Reshotho in refer- 
ence 34 has also indicated that the minimum critical Reynolds number for the 
boundary layer of an insulated plate would rise sharply with increasing Mach num- 
ber above a Mach number of 3 and indicates a need for the reexamination of the 
basic assumption of the theory of stability of the laminar boundary layer at 
high supersonic and hypersonic Mach numbers. 

Comparison of Heat-Transfer Distributions With Theory 

Distributions over smooth continuous plate .- The experimental heat-transfer 
distributions in the form of the variation of Stanton number with distance from 
the leading edge are shown in figure 12. Also, given in this figure are the cal- 
culated variations of Stanton number with distance from the leading edge for the 
laminar, transitional, and turbulent regions of the boundary layer obtained by 
the methods presented In the section entitled "Review of Heat-Transfer Equations." 
This figure clearly shows the relative values of heat transfer to be expected 
for the various types of local boundary layers with natural transition. The 
experimental Stanton number Is seen to increase by a factor of approximately 3 
from the beginning of transition to the end of transition. The calculations 
predict an increase by a factor of approximately 4; however, the leveling of the 
experimental heat-transfer values in the rearward portion of the laminar 
boundary- layer region and the choice of the virtual origin may easily account for 
the differences between the calculated and experimental values. 

Comparison of the theory with the experimental data shows that the theory 
gives a reasonably good prediction of the magnitude of heat transfer In the lam- 
inar region, generally within 15 percent of the experimental value except for the 
rearward portion of the laminar flow where the experimental heat transfer was 
found to level off prior to the rapid increase in Stanton number associated with 
the beginning of transition. 

The simple assumption of the virtual origin for the turbulent flow being 
located at the beginning of transition yields a fair prediction of the heat 
transfer In the fully developed turbulent boundary- layer region, the maximum 
deviation between theory and experiment being approximately 30 percent. Note 
particularly in figure 12(b) that the rate of decrease of turbulent Stanton num- 
ber with increasing Reynolds number follows very closely the rate of decrease 
predicted by theory. 
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The simple empirical relation for the variation of Stanton number in the 
transition region (given in eq. (15)) also leads to a reasonably good prediction 
of the experimental results. (Note that the end of transition for the 

R 0 = 5.82 X 106 and Mo = 6.0 data has been assumed to occur at 
x = 11.00 inches based on the heat-transfer-rate distribution of fig. 2 in the 
calculation of the Stanton number in the transition region.) Therefore, the 
difficulty in considering a region in which the boundary layer is transitional 
is not in determining how the heat transfer associated with the region will vary, 
but in determining where the laminar flow will end and where the fully developed 
turbulent flow will begin. 

The variation of heat transfer with local Reynolds number for the case of 
pure laminar flow is shown in figure 13 for the full test range of free-stream 
unit Reynolds number. The comparison of the experimental data with the theoret- 
ical distribution shows that theory gives a good prediction, both of the magni- 
tudes of the data and the slope of the distribution (up to the rearward region 
of laminar flow where the experimental data are observed to level off prior to 
transition) . 

Variation of Stanton number with_ Reyn olds number_ for turbulent flow.- The 
variation of the heat transfer with Reynolds number (based on the distance from 
the virtual origin for the higher Reynolds number tests of the Langley 20-inch 
Mach 6 tunnel) for the case of fully developed turbulent flow is compared with 
theory in figure 14. Also in this figure, the experimental turbulent heat- 
transfer data for the plate with roughness (for values of k equal to or greater 
than the critical height) are compared with theoretical predictions. The theory 
is seen to predict reasonably well both the magnitudes and the slopes of data for 
the plate with and without roughness. 

Inspection of figure 14 shows that the highest experimental Stanton number 
values occur near the roughness position. These peak values of Ngt decrease 
slightly as the roughness height is increased above the critical height. 

Increasing the roughness height apparently has the effect of increasing slightly 
the effective Reynolds number at any given position. However, the overall reduc- 
tion in Stanton number with increasing roughness height at any position is rather 
small since the Stanton number variation is a weak function of Reynolds number 
for turbulent flow. Therefore, it appears that roughness can be used as a 
boundary-layer trip for high-speed wind-tunnel tests designed to study heat trans- 
fer in a turbulent boundary layer without seriously affecting the heating 
distribution. 

The turbulent Stanton number variation with Reynolds number 'based on the 
distance from the roughness location for the complete test range of free-stream 
unit Reynolds number at a Mach number of 6.0 and roughness of critical height or 
greater is given in figure 15 . As explained previously, for the low Reynolds 
number tests, fully developed turbulent flow was not obtained so that the length 
of transition l could not be determined. The assumptions leading to the two 
theoretical distributions of Stanton number in figure 15 are discussed in the 
section entitled "Review of Heat-Transf er Equations." Examination of figure 15 
shows that the assumption of the virtual origin to be located at the roughness 
position gives the best agreement with the experimental data over the complete 
test Reynolds number range. 
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CONCLUSIONS 


An investigation has been conducted in the Langley 20- inch Mach 6 tunnel 
and in the Mach number 6.2 blowdown tunnel of the Langley Research Center to 
determine the effects of controlled surface roughness on boundary- layer transi- 
tion and heat transfer for local free-stream Mach numbers of 6.0 and ^-. 8 . Anal- 
ysis of the experimental results and comparison with theory and previous results 
from the literature have led to the following conclusions: 

1. A compilation of the current data and the data available from the litera- 
ture indicates that the Reynolds number required to bring about natural transi- 
tion increases markedly with increasing Mach number above a Mach number of 
approximately 3 . 7 . 

2. For a Mach number of 6.0, surface roughness that is less than the calcu- 
lated boundary- layer velocity thickness at the roughness position can under cer- 
tain conditions delay transition. 

3- The critical roughness Reynolds number required to move the beginning of 
turbulent flow to the region of the roughness is greater for higher supersonic 
and hypersonic Mach numbers than has been found at lower supersonic speeds. 

4. The high values of the required critical roughness Reynolds numbers for 
the Mach numbers of this investigation led to required roughness heights of 
approximately twice the calculated boundary- layer thickness at the roughness 
position. The effect of the large required roughness on the heat-transfer dis- 
tribution downstream of the roughness has been shown to be small. 

5- Theoretical calculations of the heat-transfer distributions agree reason- 
ably well with the experimental data obtained with the smooth flat plate. Also, 
theoretical calculations of the turbulent heat-transfer data on the plate with 
roughness of critical height or greater agreed reasonably well when it is assumed 
that fully developed turbulent flow begins at the roughness position. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va.^ October 9* 1963* 
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Figure 2.- Heating-rate distribution on the continuous flat plate. 
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Figure 3.- Concluded. 


25 




>hook from leading edge 


White band indicates 
laminar flow 











X, in. 


(c) Interchangeable leading edges; R 0 ~ 4 


X 10 6 . 


Figure 6.- Continued. 


i ■ n un ii hm 1 1 niiii min i i 1 1 i ■■ 1 1 m i n an i ■ i ■ 1 1 1 


O □ 

□ o 




(d) Interchangeable leading edges; Rq ~ 2.6 X 10 
Figure 6.- Continued. 



x, in. 


(e) Interchangeable leading edges; R 0 » 1.2 x 10°. 
Figure 6.- Continued. 
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Figure 6.- Concluded. 




Figure 7.- Heating- rate distributions on plate 2 for various size spheres. Mq = 4 . 8 . 
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Figure 15.- Comparison of experimental turbulent heat-transfer distribution (with roughness) with distribution predicted by- 

theory for complete test Reynolds number range. 
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